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FOREWORD 


This document represents the results of a developmental study by Lockheed 
ETectronics Co. - Systems and Services Division (SSD) to produce a computer 
program, which will provide energy analyses of energy conversion systems, 
such as the Modular Integrated Utility Systems (MIUS). The program title is 
Energy Systems Optimization Program 2 (ES0P2). 

A MIUS is designed to provide communities of limited size with the required 
utility services of electric power generation, heating and cooling sources, 
potable water supply, waste water treatment, and solid waste disposal. 

This study was funded by both the Department of Housing and Urban Development 
(hud) and the Systems Evaluation Office (SEO) of the NASA/Johnson Space 
Center and was under the overall direction of HUD. The technical monitor 
for the study was Mr. Richard C, Wadle of SEO. The authors wish to thank 
Mr. Wadle for his support, critique, and timely suggestions throughout the 
duration of the study. Credit is also acknowledged to other members of SEO, 
HUD, and the National Bureau of Standards, who participated in the develop- 
ment, application, and/or evaluation of the program which developed from this 
study. 

The authors also wish to acknowledge the support of SSD personnel who contri- 
buted to the development and application of this program. These include: 

Dr. J. L. Allison, who was instrumental in the incorporation of the optimiza- 
tion feature into the program, Ms. E. C. Osborne, who developed the waste 
water treatment subroutine, Mr. E. S. Riley, who incorporated the economic 
data base and financial subroutines into the program, and Mr. W. C. Rochelle 
and Mr. D. K. Liu, who were responsible for a number of analyses, using 
various applications of the program. Finally, the authors wish to acknow- 
ledge the technical support and direction of Mr. R. D. Stallings, who was 
the porject leader for two years and personally incorporated a number of 
important features into the program. 
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SUMMARY 


This document is a complete revision of previous documents describing the 
Energy Systems Optimization Program, which is used to provide analyses of 
Modular Integrated Utility Systems (MIUS). The program has now been given 
the acronym ES0P2 rather than the previously designated ESOP. Principle 
changes in the program which are described in this document include: 

•• Modifications to the input format to allow modular inputs in specified 
blocks of data. 

a An optimization feature which enables the program to search automatically 
for the minimum value of one- parameter while varying the value of other 
parameters. 

• New program option flags 'for prime mover analyses.. 

• Inclusion of solar energy for space heating and domestic hot water. 

The document is presented in four volumes: (1) Volume- 1 - Engineering Manual, 
(2) Volume II - User's Manual, (3) Volume III - ES0P2 Program Listing and 
(4) Volume IV - Economic and Financial Analysis Program and Users Guide 
The present volume (Vol. I) contains the method of analysis, description 
of program subroutines and a program flow chart. 



1. INTRODUCTION 


The Energy Systems Optimization Program was written by Lockheed Electronics 
Company/ Systems and Services Division (SSD) under the direction of the Systems 
Evaluation Office (SEO) of NASA/Johnson Space Center (OSC). The program 
was written to analyze the operating performance and costs of energy conversion 
systems, including the Modular Integrated Utility Systems (MIUS) and to compare 
the MIUS performance and costs with those of a conventional utility system. 

The Department of Housing and Urban Development (HUD) has had the overall 
direction of the MIUS Program and has sponsored, along with the SEO, the 
development and application of the Energy Systems Optimization Program 2. 

The general objectives and applications of the MIUS Program may be found in 
Ref. 3. Basically a MIUS system is designed to help alleviate the current 
energy shortage by integrating various utility services together, generating 
on-site electricity, and using prime mover recovered heat for absorption 
air conditioning, space heating, domestic hot water, waste water treatment, 
or transferring this heat to thermal storage for later use. 

The development of the Energy Systems Optimization Program as applied to MIUS 
systems began at SSD in the Fall of 1972 with the development of two sub- 
programs; solid waste disposal calculations and heating and air-conditioning 
demand calculation. The first complete user's guide was published in May 
1973 which described the analysis of the following systems integrated 
together: (1) solid waste disposal, (2) heating and cooling loads, (3) energy 
requirements, (4) electric power generation, and (5) conventional utility 
system. Further updates to the program were developed over the next eighteen 
months and incorporated such features as: multiple executions for a number 

of buildings of different types, all -electric options, potable water require- 
ments, waste water treatment, modified heating/cooling load calculations, 
thermal storage, improved low grade heat utilization, economic data base, 
financial analysis and the solar heating of buildings. 

This document describes several new and important changes to the Energy Systems 
Optimization' Program, which will hereafter be designated ES0P2. These changes 
include the following: 
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• Modifications to the input format to allow modular inputs in specified 
blocks pertaining to: (1) program option flags, (2) weather data, 

(3) building data, (4) load profiles, (5) equipment data arranged into 
varous categories (e.g., diesel/gas engines, steam power plants, solid 
waste disposal, solar energy, secondary conversion, thermal storage, 
cooling tower, and waste water treatement), (6) fuel and electricity 
data, and (7) optimization data. 

• An optimization feature which allows the program to search automatically 
for a minimum value of a parameter as a function of other system parameters. 

• New program option flags to limit the maximum number of prime movers desired, 
the number of hours at which a prime mover may operate at over 100% of its 
rated load and the maximum overload allowed on the prime mover. This 
feature is also incorporated with the optimization analysis. 

f The capability to make an energy analysis for 8,760 hours or 12 days, in 
addition to previous options. 

The present document is presented in four volumes as described below: 

• Volume I - Engineering Manual: Description of methods and types of analysis 
along with discussion of each subroutine and applications of the program. 

A flow chart of the program is also included. 

• Volume II -'User's Manual: Description of program input data and output 

data and a sample computer listing of the input and output for several 
cases. 

• Volume III - Program Listing: Complete computer listing of main program 

plus all subroutines. 

• Volume IV - Economic and Financial Analysis: Descriptive users guide of 

the economic and financial analysis program including examples. 



2. METHOD OF ANALYSIS 


2.1 ENERGY SYSTEMS CONCEPT 


The method of analysis used with the ES0P2 Program is based on the energy con- 
version system configuration shown in Figure 1. In this figure there are 
five basic subdivisions: (1) input, (2) primary sources or "prime movers", 

{3) heat storage/conversion, (4) uses, and (5) waste. Basically, energy is 
put into the system, a portion of the energy is converted to meet system 
loads and a certain amount of heat is rejected from the system. Any of 
the blocks on this figure may be used in a MIUS system; however, only the 
shaded blocks are considered for a conventional system. 

The user of ES0P2 may use the components of Figure 1 to design his own energy 
system or to evaluate a system, which may be a modified total energy or con- 
ventional system. By definition, a "conventional system", as used in this 
document, is defined as a system which: 

• Obtains the electricity to meet the domestic, auxiliary, and compression 
A/C electric demands from a fossil -fired power plant whose plant efficiency 
is input to the program (e.g., 30%). 

• Uses an all -compression A/C system (electrically driven) with an input COP 
(e.g., 4.0) to meet the space cooling demand. 

■ Uses a fossil -fired boiler to meet the domestic hot water and space heating 
demands. 

• Has no heat recovery, thermal storage, waste water treatment, solid waste 
disposal, absorption chillers, solar collectors or on-site electric genera- 
tion. 

a 

A "Total Energy System" (TES), as used in this document, is defined as a 
system, which has the classical elements, plus "MIUS" elements. The system may 


2-1 



FIGURH 1 


Energy Conversion System Configuration 
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• Use an engine/ genera tor to furnish or supplement the grid of the electricity 
for meeting all the electric demands. 

• Use recovered heat from the engine/ genera tor to meet the domestic hot water, 
space heating and space cooling (absorption chiller) demands. 

9 Use a boiler to meet any remaining heating demands not met from recovered heat. 

• Use a compression chiller to meet any remaining space cooling demands not 
met by recovered heat. 

• Have a solid waste disposal system with recovered heat available for 
meeting heating and cooling demands. 

• May have a waste water treatment facility where sludge may be used for 
input to the solid waste disposal process. 

• May have a hot and/or cold thermal storage tank which may be used as a 
secondary source of energy to the system. 

• May have an array of solar collectors, the heat from which would be used 
for meeting space heating demands and domestic hot water demands. 

As an aid in describing the method of analysis used with ES0P2, a brief 

discussion will be given for each of the blocks in Figure 1: 

• Block 1 is the commercial electric grid. All or part of the power furnished 
to the system may come from this source. The local electric rates may be 
input to the analysis and the cost of the electricity used is generated in 
the analysis. 

• Block 2 is the fuel used by the equipment making up the system. It may 
be gasoline, natural gas, coal, fuel oil, or a liquid petroleum gas (LPG). 

• Block 3 is the residential trash that is generated in the complex being 
studied. The trash is fed into the incinerator as fuel. 
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• , Block 4 is the generator which is driven by an engine. The characteris- 
tics of the generator are input parameters for the analysis. The genera- 
tor furnishes electricity to the system to supplement or replace the grid. 

• Block 5 is an engine, which drives the generator and has fuel as an input. 
Its output includes recovered heat which will be used elsewhere in the 
system. The engine characteristics are input parameters to the analysis. 
The fuel used may be diesel, gasoline, natural gas or an LPG. 

• Block 6 is the conventional boiler to produce hot water or steam. Boiler 
parameters such as pressure, temperature and fuel are input to the analysis 
and the heat generated is output to the heating requirements. 

• Block 7 represents the solar input to the analysis. Calculations are made 
to determine the solar collector energy based on the collector characteris- 
tics, orientation, and weather data for the area in which the building is 
located. The thermal energy developed by the solar collectors may be used 
in the space heating system and domestic hot water system. The prominence 
of solar collectors is expected to be greatly increased in view of the 
present energy situation. 

• Block 8 is the solid waste disposal system, which accepts residential type 
solid waste (trash), fuel, and sludge as inputs. These items are burned 
and the outputs are heat, which can be used in some other processes in the 
TES, and the honcombustible waste. 

• Block 9 is the compression air-conditioning compressors. The compressors 
are driven by electric motors and their output is chilled water and waste 
heat. 

• Block 10 is the hot thermal storage tanks that are used in TES's. Heat 
that is not needed for immediate use is stored in the tank for use later 
on in the day. This type of storage can be used to balance the 24 hour 
load on a building by storing hot water energy during low load periods and 
utilizing it during high load hours. 
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• Block n is the cold thermal storage tanks that are used in TES's. The 
cold water is stored for future use in the space cooling demand. The same 
load balancing effects can be used here as was suggested in the hot storage 
tanks . 

e Block 12 is the absorption air conditioning chiller which receives heat 
from the incinerator, engine/generator and/or boiler. The output 
is chilled water for space cooling and thermal storage and waste heat. 

• Block 13 is the domestic electric demand. This includes all electric de- 
mands other than the heating and cooling loads of the building. It would 
include all items which introduce heat into the conditioned space such 

as lighting. These load profiles are input to the analysis, based on the 
type of building being considered in the analysis. 

• Block 14 represents the auxiliary electrical demand. All electric demands 
not previously accounted for are included in the auxiliary demand. 

9 Blocks 15 and 16 are related in that weather data and building characteris- 
tics are common. The building data includes type of construction, number 
of floors, wall, roof and glass area, orientation, and general building 
parameters. The weather data is hourly data, based on local climate. The 
respective heating and cooling loads are calculated and used in the analysis. 

• Block 17 considers the hot and cold potable water requirements based on 
the type of building and its occupancy, 

• Block 18 is' the waste water treatment, which utilizes some of the recovered 
heat to produce sludge and treated water which are used in other parts of 
the system. 

’9 Block 19. is the cooling tower which dissipates the waste heat which cannot 
be used in any processes. 

9 Block 20 is the non-combustible waste from the solid waste disposal unit. 

This waste must be removed from the system, since it cannot be used for 
any other process. 



2.2 MODELING 


Both loads and equipment are modeled in an ES0P2 analysis. Section 2.2.1 pre- 
sents a description of the loads modeling, while Section 2,2.2 gives a descrip- 
tion of the equipment modeling. 

2.2.1 LOADS 

This section describes the various types of loads which are required in an 
ES0P2 analysis. Five general types of loads are modeled for each system: 
space heating and cooling, potable water, waste water, solid waste, and elec- 
tricity. If the loads are not modeled, profiles may be input. • 

The space heating and cooling loads are calculated, based on three types of 
data; weather, building, and input profiles. The weather data may be input 
by either standard weather tapes or computer cards giving indoor and outdoor 
temperatures and enthalpies, cloud cover and other information. The building 
data includes such things as latitude, longitude, time zone, azimuth angle, 
conductivity factors and surface areas of the walls, windows, and roof. 

Both single building and multiple buildings of various types (apartments, 
hotels, offices, residences, townhouses, schools and hospitals) may be modeled. 
The input 24-hour profiles include occupancy, occupant metabolic rate and ven- 
tilation rates. 

The potable water requirements may either be computed by the model if the 
building is a residential -type building or may be input in a profile if the 
building is commercial. The model will compute the hot and cold water usage 
for the kitchen, laundry, and bath for apartments, townhouse, and single- 
family residences. 

Waste water treatment may be utilized in the analysis. The change in quality 
of effluent and sludge are calculated as a function of the treatment pro- 
cesses. Presently, eleven processes are modeled. Both generic treatment 
processes and specific manufacturers equipment may be used in the model. 
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The solid waste disposal requirements are computed for three specific pro- 
cesses: incineration, pyrolysis, and combination processes. Operating para- 
meters are input and daily total energy requirements are calculated for these 
solid waste processes. 

The electric requirements are modeled by inputting hourly profiles of the 
domestic electric load and the auxiliary electric load. The domestic elec- 
tric load includes the indoor lighting and equipment electric loads, which 
introduce heat into the conditioned spaces. The auxiliary electric load 
includes the outdoor lighting and equipment electric loads, which do not in- 
troduce heat into the conditioned, spaces. 

2.2.2 EQUIPMENT 

This section describes some of the equipment characteristics modeled in an 
ES0P2 analysis. Types of prime movers which have been modeled include 
dual-fuel engines, gas turbines, and steam turbines. Size ranges investigated 
include 315-7624 KW sizes for engines and 400-7500 KW for turbines. Boilers 
have been modeled, using gas, oil, and coal as fuels. Boiler size requirements 
have ranged from hundreds to several hundred thousand lbs. of steam per hour. 
Other types of equipment modeled in analyses have included incinerators and 
pyrolysis units, absorption and compression chillers, therma-1 storage tanks, 
heat exchangers, cooling towers, and solar collectors. 


Profiles of heat recovery and fuel consumption for prime movers are input to 
the analysis as a function of generator load from zero load up to 120% load. 
Based on the' total electric demand at a particular point in time, the computer 
program interpolates on these profiles to obtain the correct value to use in 
the analysis. 

Among the operating parameters used for the equipment are the capacities of 
the units, efficiencies and coefficients of performance, and the various con- 
straints on the system. These constraints include maximum and minimum operating 
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temperatures, .pressures, flow rates, loads, number of units, and number of 
hours for overload (over 100%), 

The number of each type of equipment depends on the maximum assumed operating 
point (e.g., 90%) of the equipment. When a new unit comes on line the loads 
may be split evenly among all equipment (as in the case of prime movers), or 
a new unit may be allowed to go from its minimum operating point to its maxi- 
mum where all other equipment of the same type operate at their maximum oper- 
ating point (in the case of boilers). 

2.3 ANALYSIS TECHNIQUES • 

Three general techniques of analysis are performed in the ES0P2 Program. These 
techniques include : (1) energy analysis, (2) economic analysis, and (3) optimi- 

zation. Each of these techniques are discussed in detail in the next three 
sections. 

2.3.1 ENERGY ANALYSIS 

* 

The energy analysis performed in the ES0P2 program is divided into four 
categories: (1) space conditioning (HVAC) , (2) electric power, (3) total energy 
or MIUS Systems, and (4) conventional system. Each of these categories are 
described in the next four sections. 

2. 3. 1.1 Space Conditioning 

The space conditioning loads on a particular building are computed as the sum 
of five heat sources: (1) radiation through the windows, (2) conduction through 
the roof, walls, and windows, (3) ventilation air rate, (4) occupant metabolic 
rate, and (5) domestic electric load dissipated in the building. The detailed 
analysis of these hourly loads (including equations) may be found in the descrip- 
tion of Subroutine HEAIR. 

2. 3. 1.2 Electric Power 

The total electric power requirements are computed as the sum of (1) domestic 
electric load, (2) auxiliary electric load, (3) compression A/C electric load. 

The compression A/C electric load is calculated by an iteration method in the 
program, based on the difference between the total cooling load minus the 
cooling supplied by the absorption chiller. The description of this iteration 
method is found in the discussion of Subroutine AIRHT. 
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2. 3. 1.3 Total Energy or MI US System 


In a total energy or MIUS system, devices are used for collecting and using 
heat energy that is normally lost or wasted in a conversion process. The 
energy analysis performed in ES0P2 for such a system is designed to make 
as much use as possible of this recovered energy in meeting load demands. 

In addition the analysis permits the use of supplementary solar energy if 
the system under study is so equipped. Generally, all energy collected 
from special devices is used to meet hot water, heating, absorption air- 
conditioning. and waste treatment demands. Under certain program options 
any waste energy remaining after demands are met may be placed in thermal 
storage for use at -a later time. Otherwise it is lost. In all cases an 
hourly summary of energy input to the system as fuel, energy collected by 
various devices, energy used in meeting load demands and energy lost is 
computed. 

Recovered heat is classified as high grade or low grade depending on 
whether the temperature at which it is recovered is high enough to drive 
our absorption air-conditioner. This is not a clear cut criterion and it 
is subject to change with the state of the art in absorption chillers. How- 
ever, in ES0P2 typical sources of high grade heat are the incineration 
process and the exhaust jacket of internal combustion engines. Typical 
sources of low grade heat are solar collectors', engine oil after cooler 
and engine water jackets. Actual use of recovered heat depends on whether 
the prime mover is an internal combustion engine or a steam plant. 

2. 3. 1.3.1 Internal Combustion Engine System 

The "internal combustion engine: MIUS may, through the use of input options, 
be composed of all elements depicted in Figure 1. A minimum configuration 
win include the prime mover/ generator for meeting electric loads, one 
boiler to meet heating and hot water and absorption air-conditioning demands, 
and one of four possible air-conditioning systems. These are (!) all com- 
pressions, (2) all absorption, (3) fixed split percentage between absorption 
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and compression, and (4) floating split percentage between absorption and 
compression. If the floating split air-conditioning system is 'chosen the 
energy system may include thermal storage. 

For each hour the amount of heat recovered from the engine is determined. 

The waste disposal operating schedule is used to determine recovered heat 
available from that source. Finally, energy available from solar collectors 
is computed. Recovered heat is then allocated to meet the environmental 
control system demands and hot water demands. If more heat is available 
than is needed the remainder is either lost or is placed in thermal storage. 

If recovered heat is not sufficient to meet system demands supplementary 
heat may be taken from thermal storage if it is available. Additional 
heating requirements are satisfied by firing the boiler. Additional air- 
conditioning requirements are satisfied by either compression air-conditioning 
or by firing the boiler depending on the air-conditioning option chosen.. 

If additional compression air-conditioning is used, the additional usage of 
the engine leads to more recovered heat. This heat is considered in the 
hourly computation of the energy balance. 

2. 3. ,1.3. 2 Steam Power Plant 

Program design' limits the choice of a steam power plant to one of four 
specific systems. These systems are shown in figures 3-6 of section 3.7 
which describes the operation of subroutine STEAM. These four options 
together with a choice of three different solid waste disposal systems 
lead to 12 specific system configurations as shown in Table I. The user 
selects a specific system by specifying a mode number at porgram setup for 
both the steam power plant (options 1-4) and for solid waste disposal system. 
The modes for solid waste disposal are: (1) Incinerator, (2) Pyrolysis and 
(3) combination incinerator/pyrolysis. Methods of heat recovery and uses 
of recovered heat are depicted in figures 1-4. 
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TABLE I 

STEAM POWER PLANT CONFIGURATIONS 


mode type of system 


1 

INCINERATOR 

- 

OPERATING 

SYSTEM 

NO. ONE 

2 

PYROLYSIS 

- 

OPERATING 

SYSTEM 

NO. ONE 

3 

INCIN/PYROL 

- 

OPERATING 

SYSTEM 

NO. ONE 

4 

INCINERATOR 

- 

OPERATING 

SYSTEM 

NO. TWO 

5 

PYROLYSIS 

- 

OPERATING 

SYSTEM 

NO. TWO 

6 

INCIN/PYROL 

- 

OPERATING 

SYSTEM 

NO. TWO 

7 

INCINERATOR 

— 

OPERATING 

SYSTEM NO. THREE 

8 

PYROLYSIS 

- 

OPERATING 

SYSTEM 

NO. THREE 

9 

INCIN/PYROL 

- 

OPERATING 

SYSTEM 

NO. THREE 

10 

INCINERATOR 

- 

OPERATING 

SYSTEM 

NO. FOUR 

11 

PYROLYSIS 

- 

OPERATING 

SYSTEM 

NO. FOUR 

12 

INCIN/PYROL 

- 

OPERATING 

SYSTEM 

NO. FOUR 
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2. 3. 1.4 Conventional System 


For a conventional system the energy and fuel requirements are computed based 
on an assumed efficiency for the power plant (e.g., 305^), an all -compression 
A/C system with assumed COP (e.g., 4.0} to meet the space cooling demands and 
a boiler to meet the hot water and space heating demands. There 1s no heat 
recovery In the conventional analysis. Further discussion of the conventional 
system may be found in the description of Subroutine CONVEN- 

2.3.2 ECONOMIC ANALYSIS 

The economic analysis is performed separately, using the energy analysis 
results as input, 

I 

The discussion of the economic analyses will be divided into two general 
categories: costing analysis and financial analysis.' A simplified flow 

chart is presented in Figure 2, which shows the flow from the energy analysis 
(both design and seasonal loads) to the costing analysis and the financial 
analysis. A detailed discussion is given in volume 4, 

2. 3.2,1 Costing Analysis 

An economic data base comprising about 130 tables of various capital costs, 

0 & M costs, energy costs, debt service and administrative costs is incor- 
porated into the costing subroutine, ECOBAS. The costing of complex energy 
systems is difficult because of the factors which must be considered, 
such as: 

0 The large number of equipment components which go into making an energy 
system. 

• The regional variation in cost which may be significant in equipment, cost 
of installation or labor rates. 
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ECONOMIC/FINANCIAL ANALYSIS 



• CAPITAL COST 


• OPERATING AND MAINTENANCE COST 
® ENERGY COST 

• DEBT SERVICE 

• ADMINISTRATIVE 

PERFORM 

FINANCIAL 

ANALYSIS 

Y 

• YEARLY CASH FLOW 

• RATE OF RETURN 

• NET PRESENT VALUE 
® PAYBACK 

• YIELD 


FIGURE 2 INFORMATION FLOW THROUGH ECONOMIC/FINANCIAL ANALYSIS 
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• The significant rate of inflation in recent years makes costs estimates 
obsolete in a short time. (Data currently in ECOBAS is obsolete). 

A computerized information system was developed to aid in the costing of the 
systems. Highlights of this system include: 

(1) Developing a list of the major system components. 

(2) Developing tables of cost for labor and all major components and trans- 
ferring this data to a data base. This data base also includes city cost 
indices for both equipment and labor for 12 U. S. cities, and electricity, 
water and sewer, natural gas, oil, and coal costs for each state. 

The "equipment list", which lists the cost components of each of the major 
systems which influence the economics of a TES,is presented in Volume IV 
describing Subroutine ECOBAS. 


2. 3. 2. 2 Financial Analysis 

There are a number of financial factors which must be included in selecting 
an energy system. These include the following: 

e yearly cash flow 
0 return on investment 
» net present value 
t payback 
f yield 

Other factors to be considered are depreciation method, tax rates, system 
life and system value. The depreciation may be calculated using straight 
line, sum of years digit, or declining balance method. Further details on 
the financial analysis may be found in Volume IV in the discussion of the 
Subroutine FINANS. 
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The combination of the computerized costing and the financial analysis rou- 
tines provides a very useful tool in making investment decisions. It allows 
simulation of various system life cycles; e.g., the effect of inflation on 
the economic feasibility of the system can be studied by applying cost esca- 
lation factors to' the yearly costs. 

2.3.3 OPTIMIZATION 

The optimization capability is used to automatically drive some desired variable 
to a minimum while varying specified equipment parameters to their optimum value. 
This capability is an optional feature of the program. A single variable object 
function must be written to specify the criteria for optimization. The optimizer 
is the driving function for ES0P2 when it is utilized in the program. The sub- 
routine SEARCH defines the method used. There are 9 variables- shown in 
Table II which may presently be optimized in the program. The object 
functions to be minimized, are limited to (1) total fuel cost, (2) total 
BTU's required for heating and cooling, and (3) total purchased BTU's. 

The User's Manual (Vol . II) instructs the user On how to choose a particular 
function. 
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TABLE II OPTIMIZATION PARAMETERS 



SIZE 

na| 

mmm 

TEMP 

THICKNESS 

ORIENTATION 

PRIME MOVER A 

X 





PRIME MOVER B 

X 



* 


PRIME MOVER A 


X 




PRIME MOVER B 


X 




• 

SOLAR COLLECTORS 





X 

WALL INSULATION 




X 


BUILDING 





X 

COLD THERMAL STORAGE 

X 





HOT THERMAL STORAGE 

X 






Note: Any combination of parameters may be used in a single run with' the 
exception that size and max load for a single prime mover may not 
be used together,. 
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3. PROGRAM SUBROUTINES 


This section presents a discussion of each of the ES0P2 subroutines. A pro- 
gram flow chart showing how most of these subroutines are integrated into 
the program is given in Section 4.0. 

3.1 DRIVER MAIN 

Subroutine MAIN is the main driver of ES0P2. Subroutine MAIN calls subroutine 
LEAP and/or SEARCH depending on the optimization flag setting. 

3.2 SUBROUTINE LEAP 

In subroutine LEAP, all of the input data is read in by formatted read state- 
ments. This data is arranged in specified blocks pertaining to the following 
categories : 

• Program option flags 

• Weather data 

• Building data 

• Load profiles 

• Equipment data 

• Waste water treatment data 

• Waste disposal data 

• Fuel and electricity data 

The equipment data is further arranged into the following categories: 

• DieseT/^as engines 

• Steam power plants 

• Solid waste disposal 

• Solar energy 

• Secondary conversion systems 

• Thermal storage systems 

• Cooling tower 

• Waste water treatment 

Subroutine LEAP also prints out an echo of all of the input in three namelists 
: OUT, OUTl , and 0UT2. 


3.3 SUBROUTINE SWDP 

Subroutine SWDP calculates operating parameters and daily total energy input 
required for three types of solid waste disposal systems: incineration, pyro- 

lysis, and combination incineration/pyrolysis. Process by-products and re- 
coverable waste energy are also calculated. The processes are discussed below. 
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3.3.1 INCINERATION PROCESS 


The total daily waste (lbs) generated is a function of input values of the 
waste generation rate per person (Ibs/day), the total number of occupants 
and the sludge input from the waste water. The total moisture (lbs) and non- 
organics (lbs) are calculated from input percent (of total waste) values of 
moisture and non-organics. The total weight (lbs) of organic trash material 
(combustibles) is then calculated 

Organics (lbs) = Total Trash (lbs) - (moisture + Non-Organics) (lbs) 

The daily fuel energy requirement (BTU's) for the incineration process is com- 
posed of start-up fuel energy and supplementary fuel energy (required to sus- 
tain combustion). The total fuel energy requirement is calculated by 

Total Energy Required (BTU's) = Start-up (BTU's) + Supp. Energy Rate (BTU's/ 

HR) X Operating Time (HRS) 

where start-up and supplementary energy requirements .are input values. 

The heat generated by combustion of the trash is a function of the input 
trash heating value (BTU/lb). The total amount of heat (BTU's) generated by 
the incineration process is the sum of the fuel energy requirements, the 
trash combustion heat generation, and the energy available from sludge material 
(input). The heat recovered from the incineration process is a function of 
an input "recovery efficiency" factor. 

3.3.2 PYROLYSIS PROCESS 

The pyrolysis process is a self-sustaining solid waste disposal system (no 
start-up or supplementary fuel is required). No waste heat is recovered due 
to the nature of the process but marketable by-products are created (char, 
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tar, and combustible gas). The daily generation (lbs) of tar, char, and 
combustible gas is a function of the temperature and pressure (input) of the 
pynolysis process and the total amount of trash to be disposed of. The heating 
values of the tar, char, and gas are functions of process temperature only. 

The energy required (input) to sustain the pyrolysis process is obtained from 
combustion of the required amount of the gas by-product. The remaining por- 
tion of combustible gas, along with the tar and char by-products are assumed 
to be marketable. 

As in the incineration process, the daily operating time, size of the unit, 
and capital cost are functions of the total amount of waste to be disposed 
of (Ibs/hr). 

3.3.3 COMBINATION INCINERATION/PYROLYSIS 

The combination incineration/ pyrolysis process is actually two separate pro- 
cesses with a fixed percentage (input) of the totaT trash waste being disposed 
in each process. Waste heat is only recovered from the incineration process. 
Both incineration and pyrolysis output are printed in this mode. 

The total waste heat recovered from the incineration and combination incin- 
eration/pyrolysis processes is assumed to be available, in equal hourly in- 
crements, for utilization in meeting load demands of the MI US. The desired 
hours of utilization are 'input to this subroutine. 

3.4 SUBROUTINE HEAIR 

Subroutine HEAIR calculates the heating or cooling loads on a specific 
building with given internal and external environmental conditions for each 
hour of a 24 hour day. Input variables define the location, geometry and 
construction of the building and the internal and external environments. 

These definitions are provided by the following input variables: 

• Latitude, longitude, and time zone of the building location. 
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• On'entation 

• The areas, heat transfer coefficients, azimuth angles, tilt angles arid 
solar absorptance of all external surfaces (windows, walls and roofs). 

• Building domestic electric demand profile. 

« Building occupancy profile and the metabolic rate profile of the occupants. 

• Building ventilation rate profile. 

• Inside set-pokt dry bulb temperature profile. ^ 

$ Outdoor dry bulb temperature profile. 

■ Inside set-point enthalpy profile. 

• Outdoor air enthalpy profile. 

• Clearness of the atmosphere. 

• Amount of cloud cover.' 

• Solar reflectance of the ground surrounding the building. 

Four secondary routines are utilized by HEAIR in determining the heating and 
cooling loads. Subroutines HTLEAK, SUNl , SUN2 and SUN3 determine the heat 
transfer rate through the walls, roof and windows. SUNl calculates the sun 
rise and sun set hours for the day being used for the analysis, SUN? determines 
the direct and diffuse solar radiation fluxes for each hour of the day. These 
fluxes are used by SUNS to calculate the flux on each surface for each hour 
of the day. These fluxes are used by HTLEAK in calculating the surface tem- 
perature and heat being conducted into the building. The solar flux passing 
through the window is also determined by HTLEAK from the flux calculations 
of SUNl, SUN2 and SUNS. - The instantaneous heat gain or loss Is then converted 
into a heating or cooling Jload by HEAIR. 

The total 'heat gain or loss of the building is calculated as the sum of 
six distinct sources. 
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a) The building domestic electricity usage is totally dissipated within the 
building. Auxiliary electricity usage is not considered to be a heat load 
source. 

Q^Lec = domestic usage * 3414 BTU/KW-HR 

= BTU/HR 

b) Metabolic rate of the occupants (sensible and latent). The sensible and 
latent heat loads are calculated by empirical equations with metabolic 
rate per occupant and inside dry bulb temperature as independent variables. 
When the building is being cooled, both the sensible and latent gain are 
considered. In the heating mode only sensible is considered. 


0 

Q 


MET 

LAT 


(^SENS. ^LAT.) 

= 0.0 (Heating mode) 


No. of occupants 


c) Sensible and latent heat loads from ventilation air. 
Cooling mode: 

“ (^OUTSIDE ■ ^inside) 


Heating mode: 

^VENT ^ ^p “ (^OUTSIDE " "^INSIDE 
Where: 

w = mass flow 
H = specific enthalpy 
T = dry bulb temperature 
Cp' = specific heat of air 
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d) Conduction heat transfer from the outside environment through the walls, 
roof and windows. 


* 

Q 


COND 


^ROOF 

^ROOF 

^WALL 

'^WIND 


COND 

COND 

COND 

COND 


= UAAT 
= UAAT 


= UAAT 


U = heat transfer coefficient 
A = area 

AT = temperature differential 


e) Radiation heat transfer through the building window glass to the building 
interior. 

Qrad ' ^'window 

A = window area 

DIRSOL = direct solar radiation flux passing through the glass 
DIFSOL = diffuse solar .radiation flux passing through the glass 

_f) Provision is made for calculation of the heat gain from the use of hot 
water in the building. 

Qhw " 

Quii = Heat oain from the use of hot water in the building- 
HWF = Fraction of hot water demand that becomes a heat source 
^HWD" water demand 
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The total building heat gain or loss can be expressed as the sum of these 
five heat sources. 

%AIN " '^ELEC ‘^MET *^VENT ^COND '^RAD ^HW 

The above process is repeated for each building type specified by user input. 

The calculated total heating/ cooling demand profiles, the hot water demand 
profiles {calculated if residential -type building, input if non-residential- 
type) and the -input domestic and auxiliary electrical demand profiles are 
then summed to comprise the total 24-hour load demand profile. 

3.5 SUBROUTINE AIRHT 

Subroutine AIRHT calculates, on an hourly basis, the required energy output 
of the total energy systems (thermal storage, boiler, diesel/generator, and 
turbine/generator). The energy requirements for the steam power plant prime 
mover system are calculated internally by subroutine STEAM. 

The required energy output is determined by 

• Hot water heating loads - calculated or input 

• Domestic and auxiliary electric demands - input 

• Space heating demands - calculated 

• Air conditioning demands - calculated 

• Thermal storage system requirements - calculated 

For systems using thermal storage, AIRHT calls subroutines AIRTSS (Section 3.24) 
and HTTSS (Section 3.25). These subroutines perform part of the energy systems 
analyses. 

The domestic and auxiliary electric power demands are met by either a diesel/ 
generator or turbine/ generator system, depending on the MIUS operational mode. 
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During periods requiring heating the hot water heating demands are initially 
met (partially or fully) by low-grade waste heat energy recovered from the 
lubrication oil cooling system and water jacket cooling system (if not ebullient 
cooling) of the diesel/turbine prime mover. If there. is low-grade heat left 
over after meeting the hot water demand, it is used for the space heating 
demand. If there is enough low grade heat to satisfy both of the above 
demands, none of the high grade heat recovered from the prime mover exhaust 
and water jacket (if ebulliently cooled) and the incineration process will 
be used for these demands. If there is not enough low-grade heat to meet 
either the hot water or hot water plus space heating demands the high-grade 
recovered heat will be used. If there is not enough high-grade (plus low- 
grade) heat to satisfy these demands, the boiler will be fired to meet the 
necessary remaining demands. During non-heattng periods using the absorption 
air conditioning option, the high-grade heat is used until it is depleted or 
the total air conditioning demand has been met. If there is not enough high 
grade heat available a boiler will be fired. If the compression air conditionipg 
option has been chosen the total air conditioning demand is met by use of 
electric compressors. 

The total MIUS electric power load demand is comprised of the domestic and 
auxiliary demands (program input) and the calculated electric power required 
to drive the compress ion- type air conditioning system (when the MIUS is in 
compression A/C mode). 

The air conditioning load demands are met by an air conditioning system that 
has four modes of operation. 

1. All compression 

2. All abosrption 

3. Fixed split - compression/absorption 

4. Floating split - compression/absorption 

The "Floating Split" system operation is determined by an iterative method 
that maximizes usage of available waste heat by the absorption system and 
minimizes electrical consumption by the compression system for each hourly 
period. Thus, the ratio between compression/absorption air conditioning varies 
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hourly. The iterative method for the "Floating Ratio" system uses a technique 
involving up to 25 iterations to ensure that the sum of the absorption A/C 
load plus the compression A/C load comes within Z tons of the total air con- 
ditioning demand for the particular hour being considered. 

Subroutine AIRHT also calculates (when the low grade heat utilization option 
is selected by the program user) the operating characteristics of a specified 
cooling tower system, and the amount and utilization of 100°F waste heat 
recovered from the MIUS air-conditioning system condensers. 

The water coolant flow from the MIUS system to the cooling tower is determined 
by the heat rejection rate of the A/C condenser system, miscellaneous heat 
rejection requirements, and the cooling capacity of the cooling tower. Coolant 
water is lost through drift, evaporation, and blowdown processes. This loss 
is made up with waste water effluent flow. Requirements for the waste heat, 
before it is sent to the cooling tower, and waste water effluent flow are in- 
put to the program. Recovery, utilization, and wastage of the heat and water 
are summarized in the program output. 

The fixed-split system utilizes both compression and absorption in the ratio 
selected by the user. 

3., 6 SUBROUTINE GENRAT 


Subroutine GENRAT calculates the hourly energy requirements for prime mover/ 
generator units as a function of electrical load demands imposed. The program 
user may specify the engine/ generator set to be used, GENRAT also calculates 
the hourly rate of waste heat energy that is recovered from prime mover heat 
exchanger systems. The operating characteristics as a function of percent full 
rated load are input from subroutine LEAP. Among these operating character- 
istics are the specific fuel consumption in BTU/KW hr and the recovered waste 
heat from the engine exhaust, water jacket, and lubricating oil, all in BTU/hr'. 
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The number of specified prime mover/generator units required to meet the 
electrical load demand is determined by the criteria that a unit can only be 
loaded to an input maximum (eg. 90%) of its rated full load capacity. Thus, 
when one unit is loaded to 90% capacity, another unit is brought into operation 
and the two units operate at 45% load. The number of generators will be a 
minimum number required to meet the day's maximum auxiliary and domestic 
electric demand without exceeding a user selected maximum of generator rated 
output. Each hour's total electric demand (auxiliary, domestic and compression 
air-conditioning) will be evenly distributed among the generators. 

3.7 SUBROUTINE STEAM 


Subroutine STEAM calculates the energy requirements and operational charac- 
teristics of a HlUS-type steam power plant. Four steam power system configu- 
rations are considered and are combined with three solid waste disposal pro- 
cesses to yield twelve unique operating modes. The' three solid waste dis- 
posal processes have been presented and discussed in subroutine SWOP. Waste 
heat is available for MIUS utilization from the incinerator process and com- 
bination incinerator/ pyrolysis process. The four steam power plant configu- 
rations are discussed below. 

3.7.1 OPTION 1 

Figure 3 presents an operational schematic of the "Option 1" steam power plant 
configuration. The heating demands of the commercial /residential building are 
met by a low pressure (125 PSIG) boiler system operating at a given efficiency 
(typically 80%). The cooling and electrical demand are met by a high pressure 
(850 PSIG) boiler system which supplies, steam (850 psig, 900°F) to two turbines, 
•one of which is a straight condensing type that drives for example a 7500 KW 
generator to meet the electrical demand. The remaining turbine is also a 
straight condensing type which drives a compression- type air conditioning 
system to meet cooling demands. The operational efficiency curves of both 
turbines are considered when calculating required high pressure boiler output. 
Waste heat recovered from the waste disposal processes is assumed to be avail- 
able for utilization as part of the required high pressure boiler output. 
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3.7.2 OPTION 2 


Figure 4 presents an operational schematic of the "Option 2" steam power plant 
configuration. The electrical demands are met by a 7500 KW e.g. generator driven 
by a condensing turbine (auto-extraction type). The cooling and heating load 
demands are met by utilizing extracted steam (125 PSIG, 550° F). The cooling 
load demand is met by an absorption-type air conditioning system and the 
heating loads are met by a 100% efficient utilization as part of the total 
high pressure boiler requirements. The high pressure boiler output demand 
is calculated as a function of the electrical power demand and the steam ex- 
traction demand- The high pressure boiler fuel input demand is a function of 
the output demand and the boiler efficiency (typically 80%). 

3.7.3 OPTION 3 

Figure 5 presents an operational schematic of the "Option 3" steam power plant 
configuration. The heating and electrical demand are met by the same equip- 
ment and operational mode as described in option 2. The cooling loads are met 
by a combination compression/absorption air-conditioning system. Extracted 
steam from turbine #1 is used to drive low pressure turbine #2 (125 PS'TG in- 
let with back pressure of 12 PSIG), which in turn powers a shaft-driven com- 
pression type air conditioning system. The exhaust steam (12 PSIG) from 
turbine #2 is utilized to drive an absorption type air conditioning system. 

The air conditioning system split is determined by the user. These numbers 
may be obtained from an energy balance calculation for the operating charac- 
teristics of the system. 

3.7.4 OPTION 4 

Figure 6 presents an operational schematic of the "Option 4" steam power plant 
configuration. The only differences between options 3 and 4 is that the 
option 4 turbine #2 is driven directly by the high pressure boiler (850 psiG, 
900°F) and the air conditioning system split has a different ratio. 
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The required input to perform a steam power plant analysis consists of boiler 
efficiencies, coefficients of performance (COP) of the air conditioning system, 
energy required to obtain desired steam conditions at boiler header (Ah, enthalpy, 
BTU/lb), energy content of steam at turbine extraction points (enthalpy - BTU 
lb), and a daily profile (24 values) of heating/cooling demands (BTU/HR). 

The waste heat recovery profile is calculated by subroutine SWOP and the load 
demand profil.es are calculated by subroutine HEAIR. 

The subroutine output consists mainly of the steam power plant energy require- 
ments and the total MIUS' energy requirements (energy required for both steam 
plant and solid waste disposal operation). 
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MIUS SCHEMATIC FOR STEAM POWER PLANT PRIME MOVER 


125 PSIG 



TURBINE GEN 

#1 



FIGURE 3 STEAM POWER PLANT OPTION 1 
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MIUS SCHEMATIC FOR STEAM POWER PLANT PRIME MOVER 



FIGURE 4 STEAM POWER PLANT OPTION 2 
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MIUS SCHEMATIC FOR STEAM POWER PLANT PRIME MOVER 



FIGURE 5 STEAM POWER PLANT OPTION 3 
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MIUS SCHEMATIC FOR STEAM POWER PLANT PRIME MOVER 



ABSORPTION 


FIGURE 6 STEAM POWER PLANT -OPTION 4 








3.8 SUBROUTINE CONVEN 


Subroutine CONVEN calculates the hourly energy requirements for a conventional 
(commercial) utility power system to meet the total system load demands as 
calculated by subroutine AIRHT. The conventional utility system is composed 
of a boiler system identical to MIUS and a commercial electrical power gen- 
erating source operating at a specified input efficiency level. Only com- 
pression-type air conditioning is utilized with the COP as an input quantity 
and no waste heat energy is utilized or recovered. 

Daily, monthly, and seasonal energy requirements are calculated for compari- 
son with MIUS energy requirements. 

3.9 SUBROUTINE INTERP 

Subroutine INTERP performs linear interpolation between points of a 13 point 
two-dimensional curve. 

3.10 SUBROUTINE JNTERP 

Subroutine JNTERP performs linear interpolation between points of an 8 point 
two-dimensional curve. 

3.11 SUBROUTINE DISCOT 

Subroutine DISCOT performs single or double interpolation for continuous or 
discontinuous functions. For a detailed discussion, refer to NASA/JSC Pro- 
cedures Manual , Institutional Data Systems Division,Section 5.C.7.1 .E.2.1 . 

3.12 SUBROUTINE UNS (See 3.14) 

3:13 SUBROUTINE LAGRAN (See 3.14) 

3.14 SUBROUTINE DISSER 


Subroutines UNS, LAGRAN, and DISSER are called by subroutine DISCOT for array 
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dimension calculation and Lagrangian interpolation. Refer to DISCOT documen- 
tation for detailed discussion. 

3.15 SUBROUTINE PRINT 

Subroutine PRINT is the output print routine for subroutine HEAIR which cal- 
culates environment, domestic, and auxiliary imposed loads on the MIUS. 

3.16 SUBROUTINE OUT 

Subroutine OUT is the output print routine for subroutine SWDP which calcu- 
lates the operational parameters of the three system modes of solid waste 
disposal . 

3.17 SUBROUTINE LIST 

Subroutine LIST is the output print routine for subroutine AIRHT which cal- 
culates MIUS energy requirements. 

3.18 SUBROUTINE WATRP 

Subroutine WATRP is a water requirements subroutine in the ES0P2 Program, 
which predicts the kitchen, laundry, bath and toilet water profiles in a ty- 
pical residential-type building (i.e., single family residence, town house, 
garden-level apartment, or high-rise apartment). The subroutine determines 
the: (1) kitchen hot and cold water profiles, subdivided into those for the 
sink and dishwasher (rinse and wash water); (2) laundry hot and cold water 
profiles (rinse and wash water for the washing machine); (3) both hot and 
cold water profiles (sink and bathtub); and (4) toilet cold water profiles. 

Details of the calculations may be found in Ref. 4. 

Input consists of a residence- type flag and the number of occupants per type 
of building. The hourly, cumulative hourly, and daily requirements for the kitchet 
laundry, bath, toilet and total of the previous four are printed out. Also 
printed out is the percent of hot water to the kitchen, laundry, and bath; 
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cold water to the kitchen, laundry, bath, and toilet; and total percent of 
hot and cold water. If a non-residential type building (e.g., office 
building, regional mall, school, etc.) is being considered, subroutine 
WATPR is not used and the total heat for hot water, together with the 
hourly fractional variation, is input. 

3.19 SUBROUTINE SUNl 

Subroutine SUNl predicts the sun rise and sun set hours for the day being 
used for the analysis (see Ref. 5 and ‘6), 

3.20 SUBROUTINE SUN2 

Subroutine SUN2 calculates the direct and diffuse solar radiation fluxes 
for each hour of the day (see Ref. 5 and 6). 

3.21 SUBROUTINE SUNS 

Subroutine SUNS computes the flux on each building surface for each hour 
of the day (see Ref. 5 and 6). 

3.22 SUBROUTINE HTLEAK 

Subroutine HTLEAK predicts the surface temperature and the heat being con- 
ducted into the building. It also predicts the solar heat flux passing through 
the window from the flux calculations of SUNl, SUN2 and SUNS (see Ref. 5 and 
6 ). 
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3.23 SUBROUTINE AIRTSS 


Subroutine AIRTSS performs the air conditioning, prime mover and thermal 
storage analyses during the hours when air conditioning is required, if the 
MIUS has a thermal storage system. Total electric demand is determined and 
GENRAT is called to obtain the engine’ performance and to calculate recovered 
heat. The absorption and compression air conditioning performance is analyzed 
and the amount of cooling going into or out of cooling thermal storage is 
determined. A detailed discussion of the AIRTSS logic is presented in Ref. 7. 

3.24 SUBROUTINE HTTSS 

Subroutine HTTSS performs the thermal storage analyses when thermal storage 
is being used during the hours when the system is in a heating mode. The 
amount of energy going into or out of thermal storage is determined and the 
amount remaining in thermal storage is calculated. A detailed discussion of 
the HTTSS logic is presented in Reference 7. 

3.25 SUBROUTINE WWPROG 

Subroutine WWPROG performs the waste water treatment system analysis. The 
system water flow (equal to the potable water usage) passes through eleven 
treatment processes. The change in quality of the water as it passes through 
each process is determined. The amount of sludge recovered from the treat- 
ment is calculated and the sludge is transferred to the Waste-Disposal System 
for incineration for residential type of buildings. 
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The waste water treatment processes are logically independent in the sense 
that the order in which they are used and the number of times any one 
treatment is used is dependent only on the user's input. However, there 
may be practical limitations on the order in which the processes should be 
used. These limitations are explained fully in descriptions of the pro- 
cesses to which they apply in Reference 8. Both generic processes and 
specific manufacturer's systems are considered in the analysis. 

Subroutine WWPROG determines which processes are to be used the order in 
which they are to occur, the total head pressure involved, and the time at 
which all treatments have been processed. If there is any selection of 
equipment involved, the individual treatment subroutines print out the selec- 
tions made. If all the equipment described within a subroutine is inappro- 
priate because of input values, this information is printed out, but the 
program does not abort. It processes all remaining requests assuming no 
changes have been made in any of the waste water parameters because of the 
missing process or processes. The power requirements for each process are 
computed and a running total is calculated in the program. 

A short description of the subroutine for each of the treatment processes, 
component selection and output is given in the next sections. Each of the 
treatment subroutines computes the detention time and the change in water 
and sludge quality for that particular process. 

3:26 SUBROUTINE COMRON 

Subroutine COMPON is a component selection subroutine which calls the appro- 
priate treatment subroutine. 

3.27 SUBROUTINE RAWTRl 

Subroutine RAWTRl is the subroutine, which performs an analysis of the raw 
water pumping treatment process. The subroutine selects the optimum pump 
from a total of 25 pumps whose characteristics are built into the model. 
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3. 28 SUBROUTINE PRTRTl 


Subroutine PRTRTl is the subroutine which performs an analysis of the primary 
treatment process. The subroutine selects a clarifier tank between six. pos- 
sible types. The amount of concrete for the chosen tank is computed and the 
horsepower of the required motor is predicted. 

3.29 SUBROUTINE BIOLXl 

Subroutine BIOLXl is the subroutine which performs an analysis of the biolo- 
gical flow oxidation process. It is assumed that each bio disk is 6.5 in. in 
diameter with two sides effective in the process and that there are 17 such 
disks per running foot of unit. The minimum size is predicted, based on the 
units being able to handle either: (!) an input loading factor (gpd/ft^) or 
(2) an input detention time (hrs). 

3.30 SUBROUTINE CHMADl 

Subroutine CHMADl is the subroutine which performs an analysis of the chemi- 
cal addition/mixing process. The subroutine calculates the sizes of the units 
required, the number of each size, the required horsepower, dimensions of 
tank, amount of alum required, and amount of water to be added to the influent. 

3.31 SUBROUTINE FL0CL1 

Subroutine FLOCLl is the subroutine, which performs an analysis of the floc- 
culation treatment process. This subroutine involves only the sizing of a 
single tank for each unit, assuming a minimum detention time of 20 minutes and 
a 1:1 ratio of height to diameter. 

3.32 SUBROUTINE CHMCLl 

Subroutine CHMCLl is the subroutine which performs an analysis of the chemical 
clarification treatment process. The subroutine calculates the dynamic head 
requirement, sizing of the tank, and total amount of alum added. 
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3.33 SUBROUTINE FLOSTl 


Subroutine FLOSTl is the subroutine which performs an analysis of the flow 
stabilization treatment process. These calculations are similar to those of 
the flocculation process except that the minimum detention time is 1.75 hours. 

3.34 SUBROUTINE CRBADl 

Subroutine CRBADl is the subroutine, which performs an analysis of the carbon 
adsorption treatment process. This process assumes a minimum of two and a 
maximum of eight carbon columns, .which must be sized for each unit. It is 
assumed that a back flush equal to .75% of the 24 hour flow is subtracted from 
the effluent. 

3.35 SUBROUTINE FILTRl 

Subroutine FILTRl is the subroutine which performs an analysis of the filtra- 
tion treatment process. For this process, it is assumed that a minimum of 
two and a maximum of ten filters must be sized for each unit. A backwash of 
1.8% of the total influent is subtracted from the effluent. 

3.36 SUBROUTINE AERFLl 

Subroutine AERFLl is the subroutine which performs an analysis of the aerated 
flow equalization treatment process. The total head and pump requirements 
for this process are computed in the subroutine. 

3.37 SUBROUTINE NITRFl 

Subroutine NITRFl is the subroutine, which performs an analysis of the nitri- 
fication treatment process. The subroutine predicts the total head and pump ' 
requirements for this process. 
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3.38 SUBROUTINE OUTPUT 


Subroutine OUTPUT prints out the results of the sequential treatment of the 
waste water when all of the processes have been called. 

3.39 SUBROUTINE SESOP 

Subroutine SESOP is the executive routine for the solar energy heating and 
cooling analysis. It uses the input which is read in the main subroutine 
(LEAP) and calls the appropriate subroutines for the analysis to be performed. 


3.40 SUBROUTINE SOLENG 

Subroutine SOLENG performs the analysis of the flat plate non-tracking solar 
collectors and calculates the energy collected by the solar collectors. Sub- 
routines SUNl , SUN2, and SUN3 are called by subroutine SOLENG and are used 
to calculate the direct and diffuse solar flux on the collectors. A thermal 
balance on each of the collector surfaces is performed by SOLENG and the 
collector surface temperature is calculated. The heat transferred from the 
surface to the water is calculated for up to five collector surfaces. The 
energy generated is then used in satisfying system 

3.41 SUBROUTINE SEARCH 

Subroutine SEARCH (Reference 9) is the optimization driver. SEARCH locates 
a local minimum of a function of up to 20 variables. Each variable is sub- 
ject to an upper and lower bound. No gradients are required. The author 
of SEARCH is G. W. Westley at the Computing Technology Center, Union Carbide 
Corporation, Nuclear Division, Oak Ridge, Tennessee. 
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4. PROGRAM FLOW CHART 


READ OPTIKATtON 
OPTION FLAG 


^ IS ^ 

OPTIMIZATION 
ANALYSIS TO BE 


IS THIS ^ 
THE INITIAL 
PASS THROUGH 
the PROGRAM’ 


rREAD PROGRAM OPTION 
FLAGS, EQUIPMENT DATA, 
SOLID HASTE DISPOSAL 
DATA, WASTE WATER 
TREATMENT DATA, SOLAR 
COLLECTOR DATA, SECON- 
DARY EQUIPMENT DATA 


\PERFORHED’ 



f^EAD INITIALIZATION 
DATA fOR OPTIMIZATION 


L 

^ 


COMPAfS RESULT OF 
OBJECT FUNCTION 
WITH LAST RESULT AND 
ALTER OPTIMIZATION 
PARAMETERS 



ARE THE 
ENERGY DEMAND 
PROFILES TO BE 
. CALCULATED? ^ 


^S THE^ 
ANALYSIS TO 
DE DONE ON 
AH HOUR BY 
HOUR BASIS 
. FOR ONE ^ 
\^YEARi#^ 



WRITE TOTAL COMMUNITY 
DE^IAND PROFILES ON I'ASS 
STORAGE DEVICE FOR 4 
DAY OR 2 DAY ANALYSIS 


READ X 
' OUTDOOR > 

EHVIR0NP2NTAL 
COHOITIOHS 
HOUR BY HOUR 
FROM A 

k WEATHER / 

TAPE ^ 


WRITE OUTDOOR WEATHER 
CONDITIONS OK A HASS 
STORAGE DEVICE 


•^EAO OUTDOOR ENVIRONMENTAL 
COIIOmONS FOR EITHER A 
SEASONAL (4 DAY) ANALYSIS 
OR A DESIGII (2 DAY) 
ANALYSIS, DEPENDING ON 
OPTION FLAG INPUT. 


OPTIMIZATION^ 
TO BE DONE OH 
A BUILDING 
w PARAMETER?^ 


READ BUILDING 
PARAMETERS 


write building 
PARAMETERS ON A 
HASS STORAGE DEVICE 


COMPUTE TOTAL 
COIIMUHETY 
ENERGY LOADS 


IS THISX 
THE INITIAL 
PASS THROUGH 
^THE PROGRAM? 


has X, 
OPTIMIZATION 
BEEN COMPLETED 
. SUCCESSFULLY’ 


ARE THE ^ 
ENERGY DEMAND 
PROFILES TO BE 
s^CALCULATED’ ^ 
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IS ^ 

■^OPTIMIZATION 
ANALYSIS TO BE 
PERFORffiD? ^ 


/ HAS 

'■Optimization^ 
BEEN completed 

vSUCCESSFULLTJ . 


COMPUTE WASTE WATER 
TREATMENT VARIABLES 
AND PRINT IF FLAGGED 


WERE ANY \ 
of'twe buildings 

IN THE COMMUNITY 
FAMILY DWELLINGS^ 


^ASE THE \ 
PROFILES 
COMPUTED HOUR 
S. BY HOUR? ^ 


ENERGY LOADS PRINTED 
FOR SEASONAL OR DESIGN 
AS FLAGGED 


t MONTHLY ENERGY 
LOADS PRINTED 


IS THIS ^ 
THE INITIAL 
PASS THROUGH 
^THE PROGRAM’^ 


^ IS 'S 
OPTIMIZATION 
TO BE DONE 
ON SOLAR 
COLLECTOR 
. PARAMETERS? . 


COMPUTE HEAT 
RECOVERED 
FROM SOLAR 
COLLECTORS 


IS SPACE 

■^EATING REQUIRED^ 
^ DURING TIME 
\.C0NS1DERED>'' 


COMPUTE HIGH GRADE AND LON GRADE 
HEAT RECOVERED BY PRII'.E HOVERS 
USING A PARTICULAR GENERATOR SET 
TO SATISFY THE ELECTRIC LOAD 
REQUIREHEHTS. LIMITS KAY BE 
PLACED ON THE NUMBERS AHO SIZES 
OF PRIME HOVERS USED. 


THE HOT HATER DEMAND IS SATISFIED 
FIRST USING LOW GRADE HEAT AND 
THEN DRAHING FROM THE NIGH GRADE 
HEAT, IF THERE WAS NOT A SUfFICIEHT 
SUPPLY UF LOW GRADE HEAT. SPACE 
HEATING IS THEN SATISFIED BY 
UTILIZIKG THE REMAINING HIGH GRADE 
HEAT. IF THE REMAINING HIGH GRADE 
HEAT IS INSUFFICIENT TO MEET THE 
OEHAHO, A BOILER IS FIRED TO 
SATISFY THE SPACE HEATING DEMAND, 


^ARE'N^ 

SOLAR 
COLLECTORS 
TO BE 

,t:onsioereo?. 


COMPUTE HIGH 
GRADE HEAT FROM 
SOLID HASTE 
DISPOSAL PROGRAM 


PRINT OUTPUT 
INFORMATION 


IS THIS ^ 
THE INITIAL 
PASS THROUGH 
THE PROGRAM? 


/ IS THE ^ 
SECONDARY 
CONVERSION 
SYSTEM ALL 
STEAK OPTION 
DESIRED? > 



COMPUTE FUEL REQUIREMENTS 
FOR PRIME MOVERS AHO BOILER 


IS THE SECONDARY^ 
CONVERSION SYSTEM 
TO BE COMPRESSION 
^AIR CONOITIONING?^ 


THERMAL 
STORAGE TO BE 
s^COHSIDEREDt^ 


THE AMOUNT OF ENERGY GOING 
INTO (WASTED HIGH GRADE HEAT] 
OR OUT OF (SATISFYING THE DE- 
MAND) IS DETERMINED AND THE 
AMOUNT REMAINING IS CALCULATED 


^ IS ^ 

'OPTIHIZATIOM 
TO BE DONE 
OH SOLAR 
COLLECTOR 
Sj>ARAMETERS7^ 


CALCULATE THE ENERGY 
REQUIREMENTS AND 
OPERATIONAL CHARACTERISTICS 
OF A STEAM POWER PLAIIT 


^ IS \ 

'^PTIHIZATIO.'T' 
TO BE DONE 
ON SOLAR 
COLLECTOR 
PARAMETERS?^ 


PRINT THERMAL 

STORAGE 

INFORMATION 


^ HAS \ 
'^OPTIMIZATION 
BEEN COMPLETED 
. SUCCESSFULLY?^ 


PRINT STEAM 
POWER PLANT 
INFORMATION 


r 







■ns m SECONDARY'' 
CONVERSION SVSTEH 
TO BE ABSORPTION 
CONDITIONING?^ 


COKPOTE HIGH GRADE AND LOR GRADE 
HEAT RECOVERED BY PRIME HOVERS 
USING A PARTICULAR GENERATOR SET 
TO SATISFY THE ELECTRIC LOAD 
REQUIREMENTS. LIMITS MAY BE 
PLACED ON THE NUMBERS AND SIZES 
OF PRIHVE HOVERS USED. 


COMPUTE HIGH GRADE AND LOW GRADE 
HEAT RECOVERED BY PRIME HOVERS 
USING A PARTICULAR GENERATOR SET 
TO SATISFY THE ELEQRIC LOAD 
REQUIREMENTS. LIMITS HAY BE 
PLACED OH THE NUMBER AMD SIZES 
OF PRIME MOVERS USED. 


THE HOT HATER DEMAND IS 
SATISFIED FIRST USING LOW 
GRADE HEAT AND THEN DRAWING 
FROM THE HIGH GRADE HEAT IF 
THERE WAS HOT A SUFFICIENT 
SUPPLY OF LOU GRADF KFAT. 



IS THE SECONDARV 's 
CONVERSION SYSTEM 
TO BE A FIXED-SPLIT 
ABSORPTIOH/COHPRESSION 
- AIR CONDITIONING? - 


A PERCENTAGE OF THE COOLING 
DEMAND IS CHANGED TO COMPRESSION 
ELECTRIC DEHAND AND ADDED TO THE 
PREVIOUSLY COMPUTED BUILDING 
ELECTRIC DEMAND. 


THE HOT HATER DEMAND IS 
SATISFIED FIRST USING LOU 
GRADE MEAT AND THEN DRAWING 
FROM THE HIGH GRADE HEAT IF 
THERE WAS NOT A SUFFICIENT 
SUPPLY OF LOW GRADE HEAT. 


THE PEMAINING COOLING DEMAND 
IS SATISFIED BY USING THE 
REMAINING HIGH GRADE HEAT. 

IF THE SUPPLY IS NOT SUFFICIENT 
A BOILER IS USED TO MEET THE 
REMAINING DEMAND. 


ANY HIGH GRADE HEAT 
REiniNING IS DISCARDED 


COMPUTE FUEL REQUIRE RENTS 
FOR PRIME MOVER AND BOILEP 


IS THE SECONDARY 
CONVERSION SYSTEM 
TO BE A FLOATING-SPLIT 
ABSORPTION/COHPRESSION 
AIR CONDITIONING? ^ 


COMPUTE HIGH GRADE AND LOW GRADE 
HEAT RECOVERED BY PRIME HOVERS 
USING A PARTICULAR GENERATOR SET 
TO SATISFY THE ELECTRIC LOAD 
REQUIREUHTS. limits KAY BE 
PLACED ON THE NUMBERS AND SIZES 
OF PRIME HOVERS USED. 


THE HOT HATER DEWND IS 
SATISFIED FIRST USING LOW 
GRADE HEAT AND THEN DRAWING 
FPOP THE HIGH GRADE HEAT IF 
THERE HAS NOT A SUFFICIENT 
SUPPLY OF LOW GRADE HEAT. 


IF THE REMAINING HIGH GRADE HEAT 
IS HOT ENOUGH TO SATISFY THE 
COOLINC DEi'.PND USING ABSORPTION 
AIR CONDITIONING AH ITERATIVE 
PROCEDURE IS USED TO PRODUCE 
THE OPTIMUM COMPRESSION AIR 
CONDITIONING TO MEET THE 
COOLING DEMAND 


PROCESS COOLING 
TOWER DATA 


COMPUTE FUEL REQUIREMENTS 
FOR PRIME MOVER AND BOILER 


^ IS 'V, 
OPTIMIZATION 
ANALYSES TO BE 
PERFORMED? 


HAS 'n, 
OPTIMIZATION 
BEEN COMPLETED 
SUCCESSFULLY?.., 


THE HOT HATER DEMAND IS 
SATISFIED FIRST USING LOW 
GRADE HEAT AND THEN DRAWING 
FROM THE HIGH GRADE HEAT IF 
THERE WAS NOT A SUFFICIENT 
SUPPLY OF LOW GRADE HEAT. 


^ HAS THE ^ 
ENERGY ANALYSIS 
BEEN DONE HOUR 
BY HOUR? ^ 


ANY HIGH GRADE HEAT 
REMAINING IS DISCARDED 


ENERGY ANALYSIS 
DATA PRINTED 
MONTHLY 


ENERGY ANALYSIS DATA 
PRINTED FOR SEASONAL 
OR DESIGN AS FLAGGED 


COMPUTE FUEL REQUIREMENTS 
FOR PRIME HOVER AND BOILER 
























•5. CONCLUSION 


This document has presented an engineering description of the Energy Systems 
Optimization Program 2 (ES0P2), as used for MIUS, total energy systems, or 
conventional utility system application. A general description of the method 
of analysis was given, with a discussion of the energy systems concept, 
modeling, and analysis techniques. The details of the three analysis tech- 
niques were presented, and a description of each of the subroutines in the 
program was given. A program flow chart was presented. In summary, it is 
felt that this program is adequate to meet user needs in determining load 
and energy analyses for most types of utility or energy system desired. 

For information on the use of the program the reader is referred to 
Volume II. Volume III contains a complete program listing and Volume IV 
describes 'the operation and use of the economic data base and financial 
analysis capabilities. 
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